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THE WEAK GALERKIN METHOD FOR EIGENVALUE PROBLEMS 


HEHU XIE*, QILONG ZHAlt, AND RAN ZHANGt 


Abstract. This article is devoted to computing the eigenvalue of the Laplace eigenvalue problem 
by the weak Galerkin (WG) finite element method with emphasis on obtaining lower bounds. The 
WG method is on the use of weak functions and their weak derivatives defined as distributions. Weak 
functions and weak derivatives can be approximated by polynomials with various degrees. Different 
combination of polynomial spaces leads to different WG finite element methods, which makes WG 
methods highly flexible and efficient in practical computation. We establish the optimal-order error 
estimates for the WG finite element approximation for the eigenvalue problem. Gomparing with the 
classical nonconforming finite element method which can just provide lower bound approximation 
by linear elements with only the second order convergence, the WG methods can naturally provide 
lower bound approximation with a high order convergence (larger than 2). Some numerical results 
are also presented to demonstrate the efficiency of our theoretical results. 


Key words, weak Galerkin finite element methods, eigenvalue problem, lower bound, error 
estimate, finite element method. 
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1. Introduction. The study of eigenvalues and eigenfunctions of partial differ¬ 
ential operators both in theoretical and approximation grounds is very important in 
many helds of sciences, such as quantum mechanics, fluid mechanics, stochastic pro¬ 
cess, structural mechanics, etc. Thus, a fundamental work is to hnd the eigenvalues 
and corresponding eigenfunctions of partial differential operators. 

In this paper, we consider the following model problem: Find (A, u) such that 

{ —Alt = Am, in SI, 

M = 0, on dVl, 

Ifi = 1 , 

where is a polyhedral domain in [d = 2, 3). For simplicity, we are only concerned 
with the case d = 2, while all the conclusions can be extended to c? = 3 trivially. The 
classical variational form of problem is defined as follows: Find u S Hq{Q) and 
A S K such that b{u, u) = 1 and 

(1.2) a{u,v) = Xb{u,v), € Hq{XI), 

where 


a(u, m) = (Vm, Vm) and b{u,v) = (u,v). 


*LSEC and Institute of Computational Mathematics and Scientific/Engineering Computing, 
Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100190, China 
(hhxie@lsec.cc.ac.cn) 

^Department of Mathematics, Jilin University, Changchun, China (diqll5@mails.jlu.edu.cn). 

^Department of Mathematics, Jilin University, Changchun, China (zhangran@mail.jlu.edu.cn). 
The research of Zhang was supported in part by China Natural National Science Founda- 
tion(11271157, 11371171, 11471141), and by the Program for New Century Excellent Talents in 
University of Ministry of Education of China. 


1 



It is well known that the problem (11.21) has the eigenvalue sequence [5] 

0 < Ai < A2 < ... < Aj < ... —)■ +00 
with the corresponding eigenfunction sequence 

which satisfies the property b{ui,Uj) = Sij (<5^ denotes the Kronecker function). 

The first aim of this paper is to analyze the weak Galerkin finite element method 
for the eigenvalue problem. We will give the corresponding error estimates for the 
eigenpair approximations by the weak Galerkin finite element method based on the 
standard theory from mia. 

The eigenvalue is a positive real number, and thus it is credible if we get both the 
upper and lower bounds. In fact, a simple combination of lower and upper bounds 
will present intervals to which exact eigenvalue belongs. For the Laplace eigenvalue 
problems, since the Rayleigh quotient and minimum-maximum principle, it is easily 
to obtain the upper bounds of eigenvalue by any standard conforming finite element 
methods HHS]. For the lower bounds, the computation is of high interest and gen¬ 
erally more difficult. Influenced by the minimum-maximum principle, people try to 
obtain the lower bounds with the nonconforming element methods. In fact, the lower 
bound property of eigenvalues by nonconforming elements are observed in numeri¬ 
cal aspects at the beginning [3 Ha m |25l [301133]. After that, a series of results 
make progress in this aspect, e.g., Lin and Lin [12] proved that the non-conforming 
EQ™* rectangular element approximates exact eigenvalues associated with smooth 
eigenfunctions from below in 2006. Hu, Huang and Shen [8] gets the lower bound 
of Laplace eigenvalue problem by conforming linear and bilinear elements together 
with the mass lumping method. A general kind of expansion method, which was first 
proved in its full term in [31] by the similar argument in [1], is extensively used in 
[niiniiie] and the references cited therein. Another interesting way is provided in 
[14] and the corresponding paper cited therein. 

Our work was inspired by some recent studies of lower end approximation of 
eigenvalues by finite element discretization for some elliptic partial differential oper¬ 
ators mEKHIdlllsT]. One crucial technical ingredient that is needed in the analysis 
is some lower bound of the eigenfunction discretization error by the finite element 
method. The mainly challenge aspect is to design the high order elements to present 
the lower bound scheme. So far, the lower bound finite element methods are mainly 
first order nonconforming finite element methods. It is desired to design some high 
order numerical methods to obtain the lower bound eigenvalue approximations which 
is the second aim of this paper. 

Recently a new class of finite element methods, called weak Galerkin (WG) fi¬ 
nite element methods have been developed for the partial differential equations for its 
highly flexible and robust properties. The WG method refers to a numerical scheme 
for partial differential equations in which differential operators are approximated by 
weak forms as distributions over a set of generalized functions. It has been demon¬ 
strated that the WG method is highly flexible and robust as a numerical technique 
that employs discontinuous piecewise polynomials on polygonal or polyhedral finite 
element partitions. This thought was first proposed in [28] for a model second order 
elliptic problem in 2012, and further developed in ii[iiiiniiini[iii[iiii3iMi[33 
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with other applications. In order to enforce necessary weak continuities for approx¬ 
imating functions, proper stabilizations are employed. The main advantages of the 
WG method include the finite element partition can be of polytopal type with cer¬ 
tain regular requirements, the weak finite element space is easy to construct with 
any given approximation requirement and the WG schemes can be hybridized so that 
some unknowns associated with the interior of each element can be locally eliminated, 
yielding a system of linear equations involving much less number of unknowns than 
what it appears. 

The objective of the present paper is twofold. First, we will introduce weak 
Galerkin method for solving the Laplacian eigenvalue problem, which has the optimal- 
order error estimates. Furthermore we investigate the performance of the WG meth¬ 
ods for presenting the guaranteed lower bound approximation of the eigenvalues prob¬ 
lem under the assumption that the global mesh-size is sufficiently small. To demon¬ 
strate the potential of WG finite element methods in solving eigenvalue problems, we 
will restrict ourselves to any order WG elements (even for higher order WG hnite 
element spaces) and investigate the robustness and effectiveness of this method. 

An outline of the paper is as follows. In Section 2, the necessary notations, 
definitions of weak functions and weak derivatives are introduced. The WG hnite 
element scheme of the Laplace eigenvalue problem is stated. Error estimates for the 
boundary value problem are presented in Section 3. In Section 4, we establish the error 
estimates for the WG hnite element approximation for the eigenvalue problem. Section 
5 is devoted to presenting any higher order accuracy lower bound approximation of 
the eigenvalues. Some numerical results are presented in Section 6 to demonstrate 
the efficiency of our theoretical results and some concluding remarks are given in the 
last section. 

2. The weak Galerkin scheme for the eigenvalue problem. 

2.1. Preliminaries and notations. First, we present some notation which will 

be used in this paper. We denote (•, ■)m,u and || • \\m,u: the inner-product and the norm 
on iL™(w). If the region uj is an edge or boundary of some element, we use (•, 
instead of (•, We shall drop the subscript when to = 0 or w = 12. In this paper, 

Pr(oj) denotes the space of polynomials on w with degree no more than r. Throughout 
this paper, C denotes a generic positive constant which is independent of the mesh 
size. 

Let 7/t be a partition of the domain 12, and the elements in Tu are polygons 
satisfying the regular assumptions specified in [^. Denote by £h the edges in 7/t, 
and by the interior edges £h\d£t. For each element T G Th, hr represents the 
diameter of T, and h — maxT’gTy hx denotes the mesh size. 

2.2. A weak Galerkin scheme. Now we introduce a weak Galerkin scheme 
solving the problem (HID. For a given integer fc > 1, define the Weak Galerkin (WG) 
finite element space 

Vh = {v = [vQ,Vb) : uo|t G Pk{T),Vb\e G Pfc-i(e),Vr GTh,e G £h, and i;;, = 0 on dfl}. 

For each weak function v G Vh, we can define its weak gradient VwV by distribution 
element-wisely as follows. 
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Definition 2.1. For each v S Vh, the unique polynomial in [Pk-i{T)Y 

satisfying 

(2.1) (Vu,T,q)T =-(iio, V • q)T + • n)aT, Vq e [Pfe_i(T)]^, 

where n denotes the outward unit normal vector. 

For the aim of analysis, some projection operators are also employed in this 
paper. Let Qq denote the projection from Lf{T) onto Pk{T), Qt denote the 
projection from L^(e) onto Pfc-i(e), and Qh denote the projection from [L^(T)]^ 
onto [Pfe_i(T)]^. Combining Qq and Qfc together, we can define Qh = {Qo, Qb}, which 
is a projection from Hq{Q) onto I 4 . 

Now we define three bilinear forms on Vh that for any v,w G Vh, 

s{v,w)= ^ hf^'^'^{QbVo - Vb,QbWo - Wb)dT, 

T(^Th 

Qwiv, w) = {VwV, VwW) + s{v, w), 
b^{v,w) = (vo,wo), 

where 0 < e < 1 is a small constant parameter to be selected. With these preparations 
we can give the following weak Galerkin algorithm. 

Weak Galerkin Algorithm 1. Find Uh G Vh, Xh G M. such that b^iuh, Uh) = ^ 

and 

( 2 . 2 ) awiuh,v) = Xhbn,iuh,v), Mv GVh- 


3. Error estimates for the boundary value problem. In order to analyze 
the error of the eigenvalue problem by the weak Galerkin method, we need some 
estimates for the boundary value problem. The main idea is similar to m but some 
modifications. 


3.1. A weak Galerkin method for the Poisson equation. In this section, 
we consider the weak Galerkin method for the following Poisson equation 

(3.1) ( % 

[ w = 0, on ail. 

The corresponding weak Galerkin scheme is to find Uh G Vh such that 

(3.2) an,iuh,v) = {f,VQ), yvGVh. 


Define a semi-norm on Vh as follows 

= au,(w,w), Vue 14. 

We claim that ||| • ||| is indeed a norm on I 4 . In order to check the positive property, 
suppose |||u||| = 0 . Then we have Vu,u = 0 in T and QbVo = Vb on dT for all T G Th- 
It follows that 


(Vuo, Vuo)t = -(to, V • Vuo)t + (tq, Vuq ■ n)aT 
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= -('^ 0 , V ■ Vvo)t + {Vb, Vvo ■ n)ffT + {QbVo - Vb, Vvq ■ njgr 
= (Vvo,Vnjv)T = 0, 

so that Vq is piecewise constant and Vb = QbVo = '^o on dT. Notice that z;f, = 0 on 
dfi, we can obtain that v = 0. For the analysis, we also define another norm on as 

+ X! ^T^lIQbVo - VbWlr- 

T&Th 

Furthermore, it is easy to check that the weak Galerkin scheme (13.21) is symmetric 
and positive definite, which has a unique solution. 

The following commutative property plays an essential role in the forthcoming 
proof, which shows that the weak gradient operator is an approximation of the classical 
gradient operator. 

Lemma 3.1. ([27]) For any element T gTh, the following commutative property 
holds true, 

(3.3) Vn,iQhip) = e H\T). 


Proof. From the definition of the weak gradient (ED) and the integration by parts, 
we have that for any q € [Pk-i(T)]'^ 

{"^wiQhp), q)T = -{Qop, V • q)T + {Qbp, q • n)^ 

= -((/?, V-q)r + (<p,q-n)OT 
= (V(/3,q)T = (Qft(V(/3),q)T, 


which completes the proof. □ 

3.2. Error equation. Suppose u is the solution of (13.1|) . and Uh is the numerical 
solution of (13.21) . Denote by Ch the error that 

eh = QhU - Uh = {Qqu - Mo, QbU - Ub}. 

Then Ch should satisfy the following equation. 

Lemma 3.2. Let et be the error of the weak Galerkin scheme Then, for 

any v £Vh, we have 

(3.4) anj{eh,v) = i{u,v) + s{QhU,v), 

where 


£{u, v) = ((Vm - Qh'^u) ■n,vo- Vb)aT- 
T&Th 


Proof. From the definition of the weak gradient (ED) and the commutative prop¬ 
erty (E31) , we can obtain on each element T G Tk that 
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= -('^0, V • QftVu)r + {Vb, QhC^u) ■ n)oT 
= (V'yo, Q/tVM)T - {vo - Vb, Q/i(Vw) • n)oT 
= (Vwo, Vu)t - (Q/i(Vu) • n, Wo - Vb)dT- 

Summing over all elements and it follows that 

{VwQhU, Vu,w) = (Vwo, Vu) - ^ (Q?i(Vm) • n, wq - Vb)dT 

TGTh 

= (/) ''^o) + X! “ X! ^0 - Vb)dT 

T&Th TeTh 

= {f,vo)+^{u,v). 

Notice that the numerical solution Uh satisfies (1321). Then we can derive that 
aw{eh,v) = t{u,v) + s{QhU,v), Vw e 14, 
which completes the proof. □ 

In order to estimate the right hand side terms of (13.41) , we still need some technique 
tools introduced in [29]. 

Lemma 3.3. 1[29|1 {Trace Inequality) Let Tu be a partition of the domain LI 
into polygons in 2D or polyhedra in 3D. Assume that the partition Tt satisfies the 
Assumptions Al, A2, and A3 as stated in \29f . Let p > 1 be any real number. Then, 
there exists a constant C such that for any T G Th and edge/face e G dT, we have 

(3.5) < ChfW\9\\l,^^^ + 
for any 9 G W^'P{T). 

Lemma 3.4. (|29jl {Inverse Inequality) Let Du be a partition of the domain LI into 
polygons or polyhedra. Assume that Th satisfies all Assumptions A1 -Aj) and p > 1 be 
any real number. Then, there exists a constant C{k) such that 

(3.6) W^tWt.p < C{k)hf^ Mt,p, VT g % 

for any piecewise polynomial ip of degree no more than k on Th- 

Lemma 3.5. ([29]) Let Th he a finite element partition of LI satisfying the shape 
regularity assumptions specified in [29] and w G Il'^~^^{Lt). Then, for 0 < s < 1 we 
have 

(3.7) lk-<3ow|||_3 < 

TgTh 

(3.8) ^ h^rl^vu - Qh{Tw)\\l, < Ch^'^\\w\\U„ 

TeTh 

where C denotes a generic constant independent of mesh size h and the functions in 
the estimates. 

Suppose w G {Lt) and let ko = min{A:,fci}. With the tools above we can 

give the estimates for £{w,v) and s{QhW,v) as follows. 
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Lemma 3.6. For each element T G 7h, we have 

||Vt>o||T < C'(||V™z;||t + hj,^\\QbVo - 'y&||aT), G Vh. 
Furthermore, there is 

hr^ Iko — VbWdT < C'(||Vwo||t + hj,'^ WQbVo — VbWar), 

< CdlVuji^IlT + hrp^WQbVo - 'ybllar)- 


Proof. From the trace inequality (13.Sp and the definition of the weak gradient 
operator m, we have following inequalities for any v GVh, 

(Vno, Vvo)t = (Vwo, Vu,n)T + {QbVo - Vb, Vno • n)aT, 

< ||Vno||T||Vu,n||T + hj.'^ IIQh^^O — Vb\\dThf^\\Vvo\\aT 

< IIVuoIItIIVui^IIt + llVuollThj.^ IIQh^'o — 'Vb\\aT^ 

which implies that 

IIVuoIIt < C(||V^n||T + hpWQbVo - ^^bllar). 

Applying the Poincare inequality, we can obtain 

_1 _1 _1 
hj,^\\vo - r’bllar < h^.^ ||no - QbVoWar + hrp^\\QbVo - VbWar 

< C(||Vz;o||t + hrj.'^ lIQbUo — Ubllgr) 

< C(||Vi„n||T + hrp^WQbVo - VbWdT), 

which completes the proof. □ 

Lemma 3.7. For any v G Vh and w G the following estimates hold 

true, 

\s{QhW,v)\<Ch^'>+i\\w\\ 

III ^ III ’ 

\£{w,v)\ < C'h'=“"i||w||fco+i|||n|||, 

where ko = min{A:,/ci}. 

Proof. From the Cauchy-Schwarz inequality and Lemma 13.61 we can obtain 


\s{QhW,v)\ = 


hrp^~^'^{QbQow — QbWjQbVo — Vb)aT 


T^Th 


< C I 'y ^ hrp ~ "^^llar | ( ^ ^ ^^IIQfc^o ~ '^'bllar 

XTeTh / VreTh / 

<Ch'=»+§||n;|U„+i|lklll- 

Similarly, for the second term we can derive that 

|£(w, w)|= {{Vw-Qh'Vw) ■n,vo-Vb)dT 

T&Th 
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<C ^ ^ \\vo - ] 

XTeTh / \Terh / 

<Ch^'>-Hw\\k,+i\\vl 

which completes the proof. □ 

3.3. Error estimates. With the error equation (13.4|) and the estimates derived 
in Lemma 13.71 we can get the following error estimate for the weak Galerkin method. 

Theorem 3.8. Assume the exact solution of u e 77^+^“^ (fl), and Uh is the 
numerical solution of the weak Galerkin scheme i3.2\} . Denote ko = minjfc, fci}, then 
the following estimates hold true, 

(3.9) \lQhU-Uh\l<Ch'^°-^\\u\\ko+i, 

(3.10) IQhU - uhji < Ch'^»-^\\u\\ko+i- 


Proof. Taking v = Ch in (13.41) and it follows that 
|||e/i|f = i{u, Ch) + s{QhU, eh) 

<Ch’^«-^^\\w\\ho+4eh\l 

From the definition of ||| • |||i, we can easily get that when h is small, 
\lQhU-Uhh<Ch-^lQhU-Uh\l<Ch>^°-4u\\k,+i, 
which completes the proof. □ 

Using a standard dual argument, which is similar to the technique applied in m, 
and then we can obtain the following L^{D) error estimate. 

Theorem 3.9. Assume the exact solution of u G (fl), and Uh is the 

numerical solution of the weak Galerkin scheme Lt.A) . In addition, assume the dual 
problem has {It)-regularity. Denote ko = minjfc,/ci}, then the following estimate 
holds true 

(3.11) \\Qou-uo\\<Ch'^«+^-^\\u\\k„+i. 


4. Error estimates for the eigenvalue problem. In this section, we turn 
back to the approximation of the eigenvalue problem dm. Denote Vo = iLo(D), and 
define the sum space V = Vq + 14. Now we introduce the following semi-norm on V 
that 

II^IIe “ ^ (\\^Wo\\'^ h4\\QbWo — WbWgT^. 

Ten 

We claim that || • ||v indeed defines a norm on V. For any w £ Vo, A Wb is defined 
in the sense of trace, we shall show that Hu'lly is equivalent to |w|i, which defines a 
norm on Vq. 


Lemma 4.1. || ■ ||v^ is equivalent to | • |i on Vo- 
Proof. It is obvious that for any w GVq, 

|w|i < ||w||y. 

Then we only need to show that 

^ h^^\\Q,w-w\\lT<C\\Vw\\^. 

TdTh 

To this end, denote by Qc the projection onto Po{T), and it follows the trace 
inequality (13.51) and the Poincare’s inequality that 

hf-^WQbW- w\\Ij-< Y hf,^\\w - QcwWIj- 
T^Th TsTh 

< c Y, (^T^ll^ “ QcW||t + ||v(w - QcW)\\t) 

T&Th 

<c\\s/wr, 

which completes the proof. □ 

As to the space Vh, we have the following equivalence lemma. 

Lemma 4.2. m) There exists two constants Ci and C 2 such that for any 
w G Vh, we have 

(4.1) C'llkllv < llkllli < C2\\w\\v, 
i.e. II ■ ||y and ||| ■ |||i are equivalent on Vh- 

Proof. In order to prove the equivalence, we just need to verify the following 
inequalities that for any w GVh, 

(4.2) \\ywo\\T<C{\\V^w\\T + h~^\\Qi,wo-Wb\\aT), 

(4.3) ||Vu,w||t < C(||Vwo||t + h-j-^WQbWo - WbUdr)- 

The inequality (14.21) has been proved in Lemma 13.61 For handling the inequality (14.31) . 
we use the definition of the weak gradient to get that 

(Vu,w, Vu,ry)T = i^wW, Vwq) - {wq - Wb, V-wW ■ n)oT 

< C(||Vu,w||t||Vwo||t - \\yww\\Thj-^\\wQ - WbHar)- 

Then we can derive from Lemma 13.61 that 

IIVujwIIt < CdlVrcoIlT + hj-^\\wo - Wb\\dT) 

< C{\\Vwo\\t + hj-^ WQbWo -WbWdr), 
which completes the proof. □ 

Now we define two operators K and Kh as follows 

K : L‘^{n) ^ Vo satisfying a{Kf,v) = if, v), Mv gVo, 
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Kh : -)> Vh satisfying a^{Khf,v) = {f,Vo), € Vt- 

The following lemmas show that the finite element space 14 and the discrete 
solution operator Kh are approximations of V and K. 

Lemma 4.3. Suppose w GVq n (12), then we have 


w - Qhw\\v < Ch^^\\w\\k^+i. 


Proof. From the trace inequality (13.51) and Lemma 13.51 we have 
WQhW - w\\v 




< C 



which completes the proof. □ 

As we know, we can extend the operators K and Kh to the operators from Lf{Vl) 
to V which will not change the non-zero spectrums of the operators K and Kh- 

Lemma 4.4. The operators K and Kh have the following estimate 


lim \\Kh-K\\v = 0, 

h—^0 


where || • ||y denote the operator norm from V to V. 

Proof Since F is a Hilbert space, it is equivalent to verify that 

lim sup \\Kf - KhfWv = 0. 
ll/llv=l 

For any f GV with ||/||y = 1, suppose u = Kf and Uh = Khf. From the error 
estimate (13.101) and the regularity of the Poisson’s equation, we have 


IQhU - Uhli < Ch^-^\uh < Ch^~^\f\\ < Ch^-^WfWv. 


Then the equivalence Lemma 14.21 implies that 

WQhu - Uh\\v < 

Moreover, by letting fci = 1 in Lemma |4.3l we can obtain that 

\\u-Qhu\\v <Ch\\u\\ 2 <Ch\\f\\v. 
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It follows the triangle inequality that 

\\Kf - KhfWv = h - UhWv < Ch^-^WfWv. 

Notice that 0 < e < 1, which completes the proof. □ 

Lemma 4.5. The operator Kh '.V is compact. 

Proof. Denote the restriction of Kh on Vh. Since Vh is finite dimensional, Kh 
is compact. Notice that {Qof,vo) = {f,vo), so Kh = KhQh- In order to prove that 
Kh is compact, we just need to verify that Qh is bounded. 

For any w G 14, QhW = w. For w G Vq, we can conclude from Lemma 14.31 that 
\\Qhw\\v < WQhW - w||y + ||w||y < Cllwlli + llwllv < C'llwllv', 
which completes the proof. □ 

Now we review some notations in the spectral approximation theory. We denote 
by a{K) the spectrum of K, and by p{K) the resolvent set. Rz{K) = [z — K)~^ 
represents the resolvent operator. Let /i be a nonzero eigenvalue of T with algebraic 
multiplicities m. Let F be a circle in the complex plane centered at p, which lies in 
p{K) and encloses no other points of o-{K). The corresponding spectral projection is 

E = E{p) = ^J^Rz{K)dz. 

R{E) represents the range of E, which is the space of generalized eigenvectors. 

For a Banach space X and its closed subspaces M and N, define the distances as 
follows that 


dist(x,fV) = inf ||a: —j/||, 6{M,N) = 

y&N 


5{M, N) = max{^(M, N),5{N, M)}. 


sup dist(x,fV), 

xeM, 

||x|| = l 


Lemma 4.6. Suppose the eigenvectors of K have ki-regularity, i.e. w G 
for any w G R{E). Denote ko = min{fci,fc}, then the following estimate holds true, 


\\K-Kh\RiE)\\v<Ch'^°-^. 


Proof. Suppose w G R{E) with ||w||v = 1. Similar to the proof of Lemma lL4l 
from Theorem 13.81 we have 

WQhKw - Khw\\v < Ch^°~mTw\\ko+i. 

From Lemma 14.31 we can obtain 

\\Kw - QhKw\\v < Ch’"°\\Kw\\ko-hi, 


which implies 

\\Kw - KhwWv < Ch>^°-^\Kw\\ko+i < Ch>^°-^\w\\ko+i. 
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Since R{E) is finite dimensional, there is a uniform upper bound for ||ui||fcp+i, where 
w € R{E) with ||ru||v = 1, which completes the proof. □ 

For the symmetry of and K and Kh are self-adjoint. In addition, 

if we change the || • ||y norm to L'^(fl) norm, all the conclusions in this section can 
be interpreted trivially. Then from the theory in Bill], we can derive the following 
estimates. 

Theorem 4.7. Suppose Xj^h is the j-th eigenvalue of and Uj^h is the corre¬ 
sponding eigenvector. There exist an exact eigenvalue Xj and the corresponding exact 
eigenfunction Uj such that the following error estimates hold 

(4.4) \X, - X,,f,\ < Ch^>^o-^^\\u,U,+i, 

(4.5) \\uj - Uj^hWv < Ch^°~^\\uj\\kQ+i, 

(4.6) 11% - %>|| < 

when Uj G (fl) and ko = min{fci, k}. 

5. Lower bounds. In this section, we shall demonstrate that the approximate 
eigenvalue Xu generated by (12.21) is an asymptotic lower bound of A. About the topic 
of lower bound of the eigenvalues, please refer to [niTiiiiiin] and the references cited 
therein. In this section, the parameter s is required to be positive, i.e. 0 < e < I. 

Lemma 5.1. Suppose {X,u) is the solution of HI.2)) . and {Xh,Uh) is the solution 
of \2.2\} . Then we have the following expansion that for any v GVh, 

X- Xh = llVn - + s{uh - v,Uh - v) - A/i||mo - - A/t(||uo|P - IkoD 

-|-2(Vm - Vu,w, VwUh) - s(u, v). 


Proof. First, we have the following formulas 
a{u,u) = llVwf = A||wf, 

^wiT’ht’^h) — ||^u;^/i|| F slUfijUfi) — A^||no|| ^ 

M = IM = 1. 

Mimicking the expansion in [I] and we have the following expansion 

(Vu - VyjUh, Vu - X/wUh) 

= (Vw, Vm) + [VyjUh, ^wUh) - 2(V'U, VwUh) 

= X-GXh- 2(Vu, Vy,Uh) - s{uh,Uh) 

= X + Xh- 2(Vu - VyjV, VwUh) - 2(Vu,u, VwUh) - s{uh, Uh) 

= X + Xh- 2(Vu - Vu,w, Vy,Uh) - 2A/j(mo, %) + 2s{uh, v) - s{uh, Uh) 

= X-GXh- 2{\/u - VnjV, VwUh) + Xh{uo - vq, uq - no) - Xhiuo, uq) - Xhivo, no) 
-\-2s{uh,v) - s(uh,Uh) 

= X- Xu- 2(Vu - Vu,n, Vi,Uh) + A/,(no - % - %) + i^hiiuo, %) - (%, %)) 

-\-2s{uh,v) - s{uh,Uh). 

Rearranging the above formula and it follows that 
X-Xh = ||Vu - Vu,m|P + s{uh - v,Uh - v) - A/i||mo - vo|P - A/j(||uo|p - UnolH 
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+2(Vm - \7^v, Vy,Uh) - s{v, v), 
which completes the proof. □ 

Lemma 5.2. m) The following lower bound for the convergence rate holds for 
the exact eigenfunction u of the eigenvalue problem (HU 

||Vu-Q,,Vu|| > 


Theorem 5.3. Let Xj and Xj^h be the j-th exact eigenvalue and its corresponding 
weak Galerkin numerical approximation. Assume the corresponding eigenvector u G 
if 1+^1 (n). Denote kg = min{A:,fci}. Then if the mesh size h is small enough, there 
exists a constant C such that 


0<Xj-Xj^h<Ch^'^°-^^\\u\\ko+i- 


Proof. Take v = QhU in Lemma |5. II From the commutative property in Lemma 
13.31 there is 


X7wV = Q/jVu, 

and it follows that 

X- Xh = ||Vu - + s{uh - v,Uh - v) - A/i||uo - uo||^ - A/jdluof - Ikoll^) 

+2(Vu - X/wV, X/wUh) - s{v,v) 

= ||Vm — Q;iVu||^ + \QhU — Uh\^ — XhWQou — uqII^ — A?i(||mo||^ — IIQo^H^) 

~t“2(Vzi V u,^yjU}i^ s(^QhU, Qfluf 

Since X/wUh G [Pk-i{T)]'^, we can obtain 

(Vu - QhX/u, VwUh) = 0. 

From the error estimate (14.5p and (14.6L we have 
IQhU - Uh\f < WQhU - UhWv < 

and 

WQou-uof <Ch^’^°+^-^^\\u\\ko+l. 

Also, it follows the property of projection in Lemma 13.51 that 

WQqu - Uolp = {uo + QqU, Uq - Qou) 

= {{u - Uo) + {u- Qou), (u - Uo) - (u- Qou)) 

= \\u - uoW'^ - ||m - Qo^ll^ 

^i^Qh'^T Qh'^) — ^ ^ hrp 

TdTh 
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Table 6.1 

Convergence rates for e = 0.1 and k = 1. 


h 

1/4 

1/8 

1/16 

1/32 

1/64 

1/128 

Ai — Ai^/i 

4.6914e+0 

1.5050e+0 

4.2473e-l 

1.1518e-l 

3.0896e-2 

8.2640e-3 

order 


1.6403 

1.8251 

1.8826 

1.8984 

1.9025 

A 2 — X2,h 

2.2610e+l 

8.8734e+0 

2.7033e+0 

7.4857e-l 

2.0152e-l 

5.3843e-2 

order 


1.3494 

1.7147 

1.8525 

1.8932 

1.9041 

1 

CO 

4.4453e+l 

1.9725e+l 

6.3969e+0 

1.8123e+0 

4.9211e-l 

1.3206e-l 

order 


1.1722 

1.6246 

1.8196 

1.8808 

1.8977 

A 4 — A4^;i 

6.4193e+l 

3.1210e+l 

1.0638e+l 

3.0563e+0 

8.3024e-l 

2.2206e-l 

order 


1.0404 

1.5528 

1.7993 

1.8802 

1.9026 


— — wlloT 

T^Th 

< C/i2fco+-||ix|U„+i. 

From Lemma [5?^ we know the terms A;i||mo — Qou|P, A/i(||mo|P — ||Qom|P)), (Vm — 
QhVu, VwUh), and s{QhU, Qhu) are of higher order than || Vu — Q?iVm|P, so that 

Ch^'‘°\\u\\ko+i < \lQhU-Uhf + \\Vu-<QhVuf<Ch^'^°-^mu\\k,+i 

is the dominant term, which completes the proof. □ 

Remark 5.1. In the next section, the numerical results show that the convergence 
rates in fact tend to 2fco — e. On the other hand, the numerical eigenvalue Xh is still 
a lower hound even if e = 0. 

6. Numerical Experiments. In this section, we shall present some numerical 
results for the weak Galerkin method analyzed in the previous sections. 

6.1. Unit square domain. In the first example, we consider the problem (ED 
on the square domain 17 = (0,1)^. It has the analytic solution 

A = (m^ + n^) 7 r^, u = sin(m 7 ra::) sm{mry), 

where m, n are arbitrary integers. The first four eigenvalues are Ai = 27r^, A 2 = 
Stt^, A 3 = and A 4 = IOtt^, where A 2 and A 4 have 2 algebraic and geometric 
multiplicities. 

The uniform mesh is applied in the following examples, and h denote the mesh 
size. Different choices of the parameter e and the degree of polynomial k are presented. 
The corresponding numerical results for the first four eigenvalues are showed in Tables 
EMM From these tables, we can find the weak Galerkin method can also give the 
reasonable numerical approximations. Furthermore, the choice of e can really affect 
the convergence order which means the convergence results in Theorem 14.71 are also 
reasonable. The numerical results included in Tables 16.1116.61 shows the eigenvalue 
approximation Xj^h is a lower bound of the exact eigenvalue Xj. Tables I6.7ll6.12l 
shows the convergence behavior of the eigenfunction approximations which reveal the 
convergence results in Theorem 14.71 
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Table 6.2 

Convergence rates for e = 0.05 and k = 1. 


h 

1/4 

1/8 

1/16 

1/32 

1/64 

1/128 

Ai — \iM 

4.5174e+0 

1.3948e+0 

3.7944e-l 

9.9389e-2 

2.5770e-2 

6.6642e-3 

order 


1.6955 

1.8781 

1.9327 

1.9474 

1.9512 

A 2 — 

2.2044e+l 

8.3229e+0 

2.4379e+0 

6.5153e-l 

1.6962e-l 

4.3857e-2 

order 


1.4052 

1.7714 

1.9037 

1.9415 

1.9514 

A 3 — ^3,h 

4.3486e+l 

1.8531e+l 

5.7514e+0 

1.5674e+0 

4.1076e-l 

1.0652e-l 

order 


1.2306 

1.6880 

1.8756 

1.9320 

1.9471 

A 4 — A4_^ 

6.3226e+l 

2.9650e+l 

9.6842e+0 

2.6795e+0 

7.0366e-l 

1.8219e-l 

order 


1.0925 

1.6143 

1.8537 

1.9290 

1.9494 


Table 6.3 

Convergence rates for e = 0 and k = 1. 


h 

1/4 

1/8 

1/16 

1/32 

1/64 

1/128 

Ai — Ai^/i 

4.3486e+0 

1.2926e+0 

3.3916e-l 

8.5854e-2 

2.1531e-2 

5.3870e-3 

order 


1.7503 

1.9302 

1.9820 

1.9955 

1.9989 

A 2 — X2,h 

2.1485e+l 

7.8051e+0 

2.2003e+0 

5.6820e-l 

1.4323e-l 

3.5882e-2 

order 


1.4608 

1.8267 

1.9532 

1.9881 

1.9970 

A 3 — \3,h 

4.2518e+l 

1.7394e+l 

5.1703e+0 

1.3566e+0 

3.4342e-l 

8.6124e-2 

order 


1.2895 

1.7503 

1.9302 

1.9820 

1.9955 

A 4 — X4,h 

6.2257e+l 

2.8157e+l 

8.8226e+0 

2.3548e+0 

5.9881e-l 

1.5035e-l 

order 


1.1447 

1.6742 

1.9056 

1.9754 

1.9938 


Table 6.4 

Convergence rates for e = 0.1 and k = 2. 


h 

1/4 

1/8 

1/16 

1/32 

1/64 

1/128 

Ai — Xi,h 

2.2414e-l 

1.4066e-2 

9.2552e-4 

6.1762e-5 

4.1371e-6 

2.7718e-7 

order 


3.9941 

3.9259 

3.9055 

3.9000 

3.8997 

A 2 — X 2 M 

3.4702e+0 

1.9914e-l 

1.2641e-2 

8.3435e-4 

5.5667e-5 

3.7246e-6 

order 


4.1232 

3.9776 

3.9213 

3.9058 

3.9016 

A 3 — As^/i 

1.5809e+l 

9.7847e-l 

6.0649e-2 

3.9786e-3 

2.6525e-4 

1.7761e-5 

order 


4.0141 

4.0120 

3.9302 

3.9068 

3.9006 

A 4 — \4,h 

2.7771e+l 

2.0077e+0 

1.2222e-l 

7.9478e-3 

5.2797e-4 

3.5281e-5 

order 


3.7900 

4.0380 

3.9428 

3.9120 

3.9035 


Table 6.5 

Convergence rates for e = 0.05 and k = 2. 


h 

1/4 

1/8 

1/16 

1/32 

1/64 

1/128 

Ai — Ai^/i 

2.1062e-l 

1.2818e-2 

8.1597e-4 

5.2613e-5 

3.4037e-6 

2.2032e-7 

order 


4.0384 

3.9735 

3.9550 

3.9503 

3.9494 

A2 — A 2 ,/i 

3.2288e+0 

1.8086e-l 

1.1158e-2 

7.1301e-4 

4.5989e-5 

2.9733e-6 

order 


4.1580 

4.0188 

3.9680 

3.9546 

3.9511 

A3 — As^/i 

1.4735e+l 

8.8125e-l 

5.3250e-2 

3.3846e-3 

2.1813e-4 

1.4109e-5 

order 


4.0635 

4.0487 

3.9757 

3.9557 

3.9505 

A4 — A4./1 

2.6186e+l 

1.8074e+0 

1.0753e-l 

6.7866e-3 

4.3622e-4 

2.8175e-5 

order 


3.8568 

4.0711 

3.9859 

3.9596 

3.9525 
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Table 6.6 

Convergence rates for e = 0 and k = 2. 


h 

1/4 

1/8 

1/16 

1/32 

1/64 

1/128 

Ai — Ai^/i 

1.9798e-l 

1.1680e-2 

7.1897e-4 

4.4770e-5 

2.7957e-6 

1.7471e-7 

order 


4.0833 

4.0219 

4.0053 

4.0012 

4.0002 

A 2 — X2,h 

3.0076e+0 

1.6440e-l 

9.8518e-3 

6.0934e-4 

3.7989e-5 

2.3730e-6 

order 


4.1933 

4.0606 

4.0151 

4.0036 

4.0008 

A 3 ~ Aa^/i 

1.3731e+l 

7.9457e-l 

4.6751e-2 

2.8767e-3 

1.7911e-4 

1.1184e-5 

order 


4.1111 

4.0871 

4.0225 

4.0055 

4.0013 

A 4 — A4,;i 

2.4673e+l 

1.6299e+0 

9.4683e-2 

5.7963e-3 

3.6040e-4 

2.2497e-5 

order 


3.9201 

4.1056 

4.0299 

4.0075 

4.0018 


Table 6.7 

Convergence rates for e = 0.1 and k = 1. 


h 

1/4 

1/8 

1/16 

1/32 

1/64 

1/128 

III 

2.1026e+0 

1.1328e+0 

5.9471e-l 

3.1002e-l 

1.6128e-l 

8.3847e-2 

order 

0.0000 

0.8923 

0.9296 

0.9398 

0.9427 

0.9438 

WQhUl - 

2.4300e-l 

6.7500e-2 

1.8371e-2 

4.9711e-3 

1.3428e-3 

3.6248e-4 

order 

0.0000 

1.8480 

1.8775 

1.8858 

1.8883 

1.8893 

'iQhU2 — U2,h'i 

4.1174e+0 

2.6638e+0 

1.4918e+0 

7.9139e-l 

4.1310e-l 

2.1467e-l 

order 

0.0000 

0.6283 

0.8364 

0.9146 

0.9379 

0.9444 

\\QhU2 — U2,h\\ 

2.0401e-l 

5.0309e-2 

1.2723e-2 

3.2491e-3 

8.3336e-4 

2.1441e-4 

order 

0.0000 

2.0198 

1.9833 

1.9694 

1.9630 

1.9586 

IQ/jMs - M 3 ,hi 

4.1174e+0 

2.6638e+0 

1.4918e+0 

7.9139e-l 

4.1310e-l 

2.1467e-l 

order 

0.0000 

0.6283 

0.8364 

0.9146 

0.9379 

0.9444 

WQhUz — M3,/i 

2.0401e-l 

5.0309e-2 

1.2723e-2 

3.2491e-3 

8.3336e-4 

2.1441e-4 

order 

0.0000 

2.0198 

1.9833 

1.9694 

1.9630 

1.9586 

|||Q/i^ 4 '^4,/i||| 

8.7353e+0 

4.3751e+0 

2.3571e+0 

1.2372e+0 

6.4476e-l 

3.3534e-l 

order 

0.0000 

0.9975 

0.8923 

0.9300 

0.9402 

0.9431 

WQhUi - M4,h|| 

1.1748e+0 

2.6252e-l 

7.2940e-2 

1.9847e-2 

5.3688e-3 

1.4497e-3 

order 

0.0000 

2.1620 

1.8477 

1.8778 

1.8863 

1.8888 

lllQ/i‘^5 III 

7.6747e+0 

4.9834e+0 

2.9631e+0 

1.6029e+0 

8.4057e-l 

4.3693e-l 

order 

0.0000 

0.6230 

0.7500 

0.8864 

0.9313 

0.9440 

II Qh'^b II 

7.4542e-l 

1.1915e-l 

2.9679e-2 

7.5124e-3 

1.9132e-3 

4.8882e-4 

order 

0.0000 

2.6453 

2.0053 

1.9821 

1.9733 

1.9686 

lllQ/i‘^6 III 

6.9213e+0 

4.9834e+0 

2.9631e+0 

1.6029e+0 

8.4057e-l 

4.3693e-l 

order 

0.0000 

0.4739 

0.7500 

0.8864 

0.9313 

0.9440 

||^/i”^6 "^6, /ill 

4.8677e-l 

1.1915e-l 

2.9679e-2 

7.5124e-3 

1.9132e-3 

4.8882e-4 

order 

0.0000 

2.0304 

2.0053 

1.9821 

1.9733 

1.9686 


6.2. L shape domain. Now we consider the eigenvalue problem m on the L 
shape domain Q = (—1,1)^\(0,1)^. 

We also use the weak Galerkin method to solve this eigenvalue problem and Table 
16.13116.181 presents the corresponding numerical results for the first six eigenvalues. 
Even the analytic eigenpairs is not known, from these tables, we can hnd the numerical 
eigenvalues Xj^h increases when h decreases which shows that is a lower bound of 
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Table 6.8 

Convergence rates for e = 0.05 and k = 1. 


h 

1/4 

1/8 

1/16 

1/32 

1/64 

1/128 

WQhUl - 1 

2.0408e-b0 

1.0776e-b0 

5.5472e-l 

2.8361e-l 

1.4474e-l 

7.3831e-2 

order 


0.9213 

0.9580 

0.9678 

0.9704 

0.9712 

WQhUl - Ml, /ill 

2.2929e-l 

6.1245e-2 

1.6035e-2 

4.1758e-3 

1.0858e-3 

2.8224e-4 

order 


1.9045 

1.9333 

1.9411 

1.9432 

1.9438 

\lQhU2 — U2,h\l 

4.0525e-b0 

2.5606e-b0 

1.4016e-b0 

7.2851e-l 

3.7309e-l 

1.9031e-l 

order 


0.6623 

0.8693 

0.9441 

0.9654 

0.9711 

\\QhU2 — U2,h\\ 

2.0213e-l 

4.9226e-2 

1.2287e-2 

3.0925e-3 

7.8035e-4 

1.9713e-4 

order 


2.0378 

2.0023 

1.9903 

1.9866 

1.9850 

iQ/iM3 - M3,/i 1 

4.0525e-b0 

2.5606e-b0 

1.4016e-b0 

7.2851e-l 

3.7309e-l 

1.9031e-l 

order 


0.6623 

0.8693 

0.9441 

0.9654 

0.9711 

\\QhU3 — M3,/i 

2.0213e-l 

4.9226e-2 

1.2287e-2 

3.0925e-3 

7.8035e-4 

1.9713e-4 

order 


2.0378 

2.0023 

1.9903 

1.9866 

1.9850 


7.4014e+0 

4.1637e-b0 

2.1988e-b0 

1.1318e-b0 

5.7863e-l 

2.9528e-l 

order 


0.8299 

0.9211 

0.9581 

0.9679 

0.9705 

\\QhU4. - U4,h\\ 

8.7222e-l 

2.3834e-l 

6.3677e-2 

1.6672e-2 

4.3413e-3 

1.1288e-3 

order 


1.8717 

1.9042 

1.9333 

1.9412 

1.9433 

lllQ/i‘^5 III 

6.8886e-b0 

4.8293e-b0 

2.8003e-b0 

1.4825e-b0 

7.6263e-l 

3.8924e-l 

order 


0.5124 

0.7862 

0.9176 

0.9590 

0.9703 

\\QhU5 — M5,/i 

4.9173e-l 

1.1747e-l 

2.8947e-2 

7.2442e-3 

1.8217e-3 

4.5881e-4 

order 


2.0656 

2.0208 

1.9985 

1.9916 

1.9893 

lllQ/i‘^6 III 

6.8882e-b0 

4.8293e-b0 

2.8003e-b0 

1.4825e-b0 

7.6263e-l 

3.8924e-l 

order 


0.5123 

0.7862 

0.9176 

0.9590 

0.9703 

II Qh'^6 II 

4.9159e-l 

1.1747e-l 

2.8947e-2 

7.2442e-3 

1.8217e-3 

4.5881e-4 

order 


2.0652 

2.0208 

1.9985 

1.9916 

1.9893 


the exact eigenvalue Xj. 

7 . Concluding remarks. In this paper, we apply the weak Galerkin method 
to solve the eigenvalue problems and the corresponding convergence analysis is also 
given. Furthermore, we also analyze the lower-bound property of the weak Galerkin 
method. Gompared with the classical nonconforming finite element method which 
can provide lower bound approximation by linear element with only the second order 
convergence, the weak Galerkin method can provide lower bound approximation with 
a high order convergence (larger than 2). 

In the future, it is required to design the efficient solver for the algebraic eigenvalue 
problems derived by the weak Galerkin method. 
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Table 6.9 

Convergence rates for e = 0 and k = 1. 


h 

1/4 

1/8 

1/16 

1/32 

1/64 

1/128 

ilQ/i^i ^i,/i||| 

1.9806e+0 

1.0247e+0 

5.1693e-l 

2.5905e-l 

1.2960e-l 

6.4808e-2 

order 


0.9508 

0.9871 

0.9967 

0.9992 

0.9998 

WQhUi - ui^hW 

2.1635e-l 

5.5551e-2 

1.3986e-2 

3.5027e-3 

8.7607e-4 

2.1904e-4 

order 


1.9615 

1.9898 

1.9974 

1.9994 

1.9998 

lQhU2 - U2.hl 

3.9877e+0 

2.4601e+0 

1.3162e+0 

6.7019e-l 

3.3666e-l 

1.6852e-l 

order 


0.6968 

0.9024 

0.9737 

0.9933 

0.9983 

\\QhU2 — U2,h\\ 

2.0039e-l 

4.8309e-2 

1.1944e-2 

2.9773e-3 

7.4379e-4 

1.8591e-4 

order 


2.0524 

2.0160 

2.0042 

2.0011 

2.0003 

IQhUs - Ms,/ill 

3.9877e+0 

2.4601e+0 

1.3162e+0 

6.7019e-l 

3.3666e-l 

1.6852e-l 

order 


0.6968 

0.9024 

0.9737 

0.9933 

0.9983 

WQhUz - M3,/i 

2.0039e-l 

4.8309e-2 

1.1944e-2 

2.9773e-3 

7.4379e-4 

1.8591e-4 

order 


2.0524 

2.0160 

2.0042 

2.0011 

2.0003 

|||Q/i^4 1^4, If 

1.2992e+l 

3.9612e+0 

2.0494e+0 

1.0339e+0 

5.1810e-l 

2.5920e-l 

order 


1.7136 

0.9508 

0.9871 

0.9967 

0.9992 

WQhUi - Ui,hW 

1.9731e+0 

2.1635e-l 

5.5551e-2 

1.3986e-2 

3.5027e-3 

8.7607e-4 

order 


3.1890 

1.9615 
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1.1603e-l 
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2.0319 

2.0090 

2.0023 

2.0006 
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6.9191e-l 
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0.8223 
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0.9865 
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9.9163e-l 

1.1603e-l 

2.8373e-2 
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1.7595e-3 
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2.0006 
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Table 6.10 

Convergence rates for e = 0.1 and k = 2. 


h 

1/4 

1/8 

1/16 

1/32 

1/64 

1/128 

’^l,h III 

4.1653e-l 

1.0911e-l 

2.8410e-2 

7.3891e-3 

1.9213e-3 

4.9944e-4 

order 


1.9327 

1.9413 

1.9429 

1.9433 

1.9437 

WQhUi - ui^hW 

8.5534e-2 

1.0676e-2 

1.4001e-3 

1.8702e-4 

2.5108e-5 

3.3748e-6 

order 


3.0021 

2.9308 

2.9043 

2.8970 

2.8953 

lQhU2 - U2,hl 

1.5575e+0 

4.3076e-l 

1.1323e-l 

2.9496e-2 

7.6668e-3 

1.9916e-3 

order 


1.8543 

1.9277 

1.9406 

1.9438 

1.9447 

\\QhU2 — U2,h\\ 

3.2591e-l 

4.1406e-2 

5.2583e-3 

6.9630e-4 

9.3495e-5 

1.2600e-5 

order 


2.9766 

2.9772 

2.9168 

2.8967 

2.8914 

IQhUs - Ms.hll 

1.7381e+0 

4.4628e-l 

1.1430e-l 

2.9567e-2 

7.6714e-3 

1.9919e-3 

order 


1.9615 

1.9651 

1.9508 

1.9464 

1.9453 

WQhUz - 

3.5920e-l 

4.5291e-2 

5.7680e-3 

7.6165e-4 

1.0170e-4 

1.3623e-5 

order 


2.9875 

2.9731 

2.9209 

2.9047 

2.9003 

||iQ/i1^4 1^4,III 

3.1633e+0 

8.8821e-l 

2.2815e-l 

5.9192e-2 

1.5378e-2 

3.9962e-3 

order 


1.8324 

1.9609 

1.9465 

1.9446 

1.9441 

WQhUi - Ui^hW 

6.6471e-l 

9.6376e-2 

1.1636e-2 

1.5123e-3 

2.0146e-4 

2.7020e-5 

order 


2.7860 

3.0501 

2.9438 

2.9082 

2.8984 

|||^5/i1^5 1^5,III 

3.3748e+0 

1.2150e+0 

3.2184e-l 

8.4035e-2 

2.1843e-2 

5.6714e-3 

order 


1.4738 

1.9166 

1.9373 

1.9438 

1.9454 

WQhUh — Uz^hW 

6.5236e-l 

1.1991e-l 

1.4765e-2 

1.9227e-3 

2.5710e-4 

3.4647e-5 

order 


2.4437 

3.0217 

2.9410 

2.9027 

2.8915 

|||Q/i1^6 1^6,III 

4.2203e+0 

1.3252e+0 

3.2987e-l 

8.4571e-2 

2.1878e-2 

5.6737e-3 

order 


1.6711 

2.0063 

1.9637 

1.9507 

1.9471 

WQhUe — UQ^hW 

7.2192e-l 

1.3575e-l 

1.6679e-2 

2.1691e-3 

2.8817e-4 

3.8522e-5 

order 


2.4109 
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Table 6.11 

Convergence rates for s = 0.05 and k = 2. 


h 

1/4 

1/8 

1/16 

1/32 

1/64 

1/128 

lll,h III 

4.0028e-l 

1.0330e-l 

2.6433e-2 

6.7483e-3 

1.7218e-3 

4.3922e-4 

order 


1.9542 

1.9664 

1.9698 

1.9706 

1.9709 

WQhUi - ui^hW 

7.9910e-2 

9.6767e-3 

1.2274e-3 

1.5829e-4 

2.0503e-5 

2.6586e-6 

order 


3.0458 

2.9789 

2.9550 

2.9486 

2.9471 

lQhU2 - U2,hl 

1.4826e+0 

4.0594e-l 

1.0531e-l 

2.6977e-2 

6.8859e-3 

1.7561e-3 

order 


1.8688 

1.9466 

1.9649 

1.9700 

1.9713 

\\QhU2 — U2,h\\ 

2.9953e-l 

3.7029e-2 

4.5641e-3 

5.8360e-4 

7.5522e-5 

9.8044e-6 

order 


3.0160 

3.0203 

2.9673 

2.9500 

2.9454 

IQhUs - Ms.hll 

1.6545e+0 

4.2057e-l 

1.0631e-l 

2.7041e-2 

6.8900e-3 

1.7563e-3 

order 


1.9760 

1.9840 

1.9751 

1.9726 

1.9719 

WQhUz - 

3.3275e-l 

4.0902e-2 

5.0762e-3 

6.4953e-4 

8.3845e-5 

1.0847e-5 

order 


3.0242 

3.0104 

2.9663 

2.9536 

2.9504 

liQhlll ll4,h||| 

3.0151e+0 

8.3014e-l 

2.1137e-l 

5.3984e-2 

1.3775e-2 

3.5138e-3 

order 


1.8608 

1.9736 

1.9692 

1.9705 

1.9709 

WQhUi - Ui^hW 

6.2435e-l 

8.5941e-2 

1.0140e-2 

1.2774e-3 

1.6440e-4 

2.1281e-5 

order 


2.8609 

3.0833 

2.9888 

2.9579 

2.9496 

|||^5hll5 III 

3.2512e+0 

1.1345e+0 

2.9867e-l 

7.6888e-2 

1.9646e-2 

5.0100e-3 

order 


1.5189 

1.9255 

1.9577 

1.9685 

1.9713 

WQhUh — Uz^hW 

6.0052e-l 

1.0549e-l 

1.2715e-2 

1.6035e-3 

2.0673e-4 

2.6825e-5 

order 


2.5090 

3.0525 

2.9873 

2.9554 

2.9461 

|||Qhll6 III 

4.0657e+0 

1.2374e+0 

3.0612e-l 

7.7378e-2 

1.9678e-2 

5.0120e-3 

order 


1.7162 

2.0152 

1.9841 

1.9754 

1.9731 

WQhUe — UQ^hW 

6.8404e-l 

1.2118e-l 

1.4639e-2 

1.8523e-3 

2.3824e-4 

3.0778e-5 

order 


2.4969 

3.0493 

2.9824 

2.9588 

2.9525 
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Table 6.12 

Convergence rates for e = 0 and k = 2. 


h 

1/4 

1/8 

1/16 

1/32 

1/64 

1/128 

IQhUl - 1 

3.8476e-l 

9.7751e-2 

2.4561e-2 

6.1490e-3 

1.5379e-3 

3.8452e-4 

order 


1.9768 

1.9927 

1.9980 

1.9994 

1.9998 

WQhUl - Ml, /ill 

7.4702e-2 

8.7708e-3 

1.0752e-3 

1.3375e-4 

1.6701e-5 

2.0872e-6 

order 


3.0904 

3.0281 

3.0070 

3.0015 

3.0003 

\lQhU2 — M2,/i 

1.4123e+0 

3.8267e-l 

9.7879e-2 

2.4631e-2 

6.1684e-3 

1.5428e-3 

order 


1.8839 

1.9670 

1.9906 

1.9975 

1.9993 

\\QhU2 — U2,h\\ 

2.7564e-l 

3.3148e-2 

3.9609e-3 

4.8855e-4 

6.0870e-5 

7.6034e-6 

order 


3.0558 

3.0650 

3.0192 

3.0047 

3.0010 

IQhUs - M3,/i 1 

1.5761e+0 

3.9645e-l 

9.8807e-2 

2.4690e-2 

6.1721e-3 

1.5431e-3 

order 


1.9911 

2.0045 

2.0007 

2.0001 

2.0000 

WQhUs — M3,/i 

3.0857e-l 

3.6994e-2 

4.4723e-3 

5.5453e-4 

6.9212e-5 

8.6507e-6 

order 


3.0603 

3.0482 

3.0117 

3.0022 

3.0001 

\lQhU4 - U4,h\l 

2.8709e+0 

7.7652e-l 

1.9567e-l 

4.9133e-2 

1.2299e-2 

3.0759e-3 

order 


1.8864 

1.9886 

1.9937 

1.9981 

1.9995 

\\QhU4 - U4,h\\ 

5.8517e-l 

7.6776e-2 

8.8369e-3 

1.0776e-3 

1.3385e-4 

1.6705e-5 

order 


2.9301 

3.1190 

3.0358 

3.0091 

3.0022 

lllQ/i'^5 '^5,h III 

3.1281e+0 

1.0609e+0 

2.7714e-l 

7.0267e-2 

1.7634e-2 

4.4131e-3 

order 


1.5600 

1.9366 

1.9797 

1.9945 

1.9985 

WQhUz — M5,/i 

5.5149e-l 

9.3075e-2 

1.0962e-2 

1.3372e-3 

1.6608e-4 

2.0730e-5 

order 


2.5669 

3.0858 

3.0352 

3.0093 

3.0021 

lllQ/i'^6 III 

3.9117e+0 

1.1571e+0 

2.8405e-l 

7.0715e-2 

1.7663e-2 

4.4148e-3 

order 


1.7573 

2.0263 

2.0061 

2.0013 

2.0003 

II QhUQ II 

6.4688e-l 

1.0852e-l 

1.2881e-2 

1.5861e-3 

1.9765e-4 

2.4697e-5 

order 


2.5756 

3.0746 

3.0217 

3.0044 

3.0005 


Table 6.13 

Discrete eigenvalues for £ = 0.1 and k = 1. 


h 

1/4 

1/8 

1/16 

1/32 

1/64 

1/128 

Ai,/i 

8.0444 

9.1197 

9.4787 

9.5893 

9.6234 

9.6343 

A2,/i 

12.1745 

14.2566 

14.9345 

15.1263 

15.1782 

15.1922 

As./i 

15.0478 

18.2342 

19.3145 

19.6240 

19.7083 

19.7309 

M,h 

19.7958 

26.1541 

28.5461 

29.2554 

29.4501 

29.5024 

As./i 

20.2283 

27.6970 

30.6403 

31.5476 

31.8078 

31.8819 

Ae./i 

23.7850 

34.6725 

39.3995 

40.8900 

41.3124 

41.4292 


Table 6.14 

Discrete eigenvalues for e = 0.05 and k = 1. 


h 

1/4 

1/8 

1/16 

1/32 

1/64 

1/128 

Ai,/i 

8.0934 

9.1472 

9.4896 

9.5931 

9.6247 

9.6346 

A2,/i 

12.2877 

14.3238 

14.9615 

15.1356 

15.1813 

15.1931 

As./i 

15.2218 

18.3444 

19.3598 

19.6398 

19.7134 

19.7325 

M,h 

20.1002 

26.3819 

28.6452 

29.2905 

29.4615 

29.5060 

As./i 

20.5464 

27.9527 

30.7545 

31.5885 

31.8212 

31.8861 

Ae./i 

24.2289 

35.0750 

39.5886 

40.9586 

41.3349 

41.4363 
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Table 6.15 

Discrete eigenvalues for £ = 0 and k = 1. 


h 

1/4 

1/8 

1/16 

1/32 

1/64 

1/128 

M,h 

8.1405 

9.1725 

9.4992 

9.5963 

9.6257 

9.6349 

^2,h 

12.3971 

14.3860 

14.9856 

15.1437 

15.1838 

15.1939 

^3,h 

15.3906 

18.4466 

19.4000 

19.6534 

19.7177 

19.7338 

M,h 

20.3977 

26.5942 

28.7336 

29.3206 

29.4710 

29.5088 

^5,h 

20.8576 

28.1912 

30.8564 

31.6235 

31.8322 

31.8894 


24.6655 

35.4522 

39.7576 

41.0175 

41.3535 

41.4419 


Table 6.16 

Discrete eigenvalues for e = 0.1 and k = 2. 


h 

1/4 

1/8 

1/16 

1/32 

1/64 

1/128 

M,h 

9.5538 

9.6152 

9.6306 

9.6362 

9.6383 

9.6392 

^2,h 

15.0957 

15.1903 

15.1967 

15.1972 

15.1972 

15.1973 

^3,h 

19.5148 

19.7251 

19.7383 

19.7391 

19.7392 

19.7392 

M,h 

28.7653 

29.4755 

29.5185 

29.5213 

29.5215 

29.5215 

^5,h 

30.6553 

31.7870 

31.8860 

31.9036 

31.9091 

31.9113 


39.0485 

41.2933 

41.4490 

41.4674 

41.4719 

41.4735 


Table 6.17 

Discrete eigenvalues for e = 0.05 and k = 2. 


h 

1/4 

1/8 

1/16 

1/32 

1/64 

1/128 

Ai,h 

9.5555 

9.6155 

9.6307 

9.6362 

9.6383 

9.6392 

^2,h 

15.1012 

15.1908 

15.1968 

15.1972 

15.1972 

15.1973 

^3,h 

19.5283 

19.7264 

19.7384 

19.7392 

19.7392 

19.7392 

M,h 

28.8132 

29.4796 

29.5188 

29.5213 

29.5215 

29.5215 

^5,h 

30.7260 

31.7934 

31.8867 

31.9037 

31.9092 

31.9113 

Ae.h 

39.2028 

41.3060 

41.4501 

41.4675 

41.4719 

41.4735 


Table 6.18 

Discrete eigenvalues for £ = 0 and k = 2. 


h 

1/4 

1/8 

1/16 

1/32 

1/64 

1/128 

Ai,ft 

9.5571 

9.6157 

9.6308 

9.6362 

9.6383 

9.6392 

^2,h 

15.1064 

15.1913 

15.1968 

15.1972 

15.1972 

15.1973 

^3,h 

19.5410 

19.7275 

19.7385 

19.7392 

19.7392 

19.7392 

M,h 

28.8577 

29.4833 

29.5191 

29.5213 

29.5215 

29.5215 

^5,h 

30.7916 

31.7993 

31.8873 

31.9038 

31.9092 

31.9113 

Ae.h 

39.3450 

41.3176 

41.4512 

41.4676 

41.4719 

41.4735 
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